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Metamaterials  exhibit  numerous  novel  effects*"^  and  operate 
over  a  large  portion  of  the  electromagnetic  spectrum^*". 
Metamaterial  devices  based  on  these  effects  include  gradient- 
index  lenses**’*^,  modulators  for  terahertz  radiation^^"^^  and 
compact  waveguides^*.  The  resonant  nature  of  metamaterials 
results  in  frequency  dispersion  and  narrow  bandwidth 
operation  where  the  centre  frequency  is  fixed  by  the  geometry 
and  dimensions  of  the  elements  comprising  the  metamaterial 
composite.  The  creation  of  frequency-agile  metamaterials 
would  extend  the  spectral  range  over  which  devices  function 
and,  further,  enable  the  manufacture  of  new  devices  such  as 
dynamically  tunable  notch  filters.  Here,  we  demonstrate  such 
frequency-agile  metamaterials  operating  in  the  far-infrared  by 
incorporating  semiconductors  in  critical  regions  of  metallic 
split-ring  resonators.  For  this  first-generation  device,  external 
optical  control  results  in  tuning  of  the  metamaterial  resonance 
frequency  by  ~20%.  Our  approach  is  integrable  with  current 
semiconductor  technologies  and  can  be  implemented  in  other 
regions  of  the  electromagnetic  spectrum. 

Electromagnetic  metamaterials  are  structured  composites  with 
patterned  metallic  subwavelength  inclusions.  These  mesoscopic 
systems  are  built  from  the  bottom  up,  at  the  unit  cell  level,  to 
yield  specific  electromagnetic  properties.  Individual  components 
respond  resonantly  to  the  electric,  magnetic  or  both  components 
of  the  electromagnetic  field.  In  this  way  electromagnetic 
metamaterials  can  be  designed  to  yield  a  desired  response  at 
frequencies  from  the  microwave  through  to  the  near  visible*"'®. 
Importantly,  additional  design  flexibility  is  afforded  by  the 
judicious  incorporation  of  naturally  occurring  materials  within 
or  as  part  of  the  metamaterial  elements.  Specifically,  hybrid 
metamaterial  composites  result  when  the  properties  of  a  natural 
material,  for  example,  semiconductors,  strongly  couple  with  the 
resonance  of  a  metamaterial  element.  The  resulting  hybrid 
metamaterials  will  still  exhibit  ‘passive’  properties  (such  as 
a  negative  electric  response,  negative  index  or  gradient  index),  as 
determined  by  the  patterning  of  the  metamaterial  elements. 
However,  the  aforementioned  coupling  engenders  control  of  the 
passive  metamaterial  response  by  means  of  an  external  stimulus 
of  the  natural  material  response  (such  as  photoconductivity, 
nonlinearity,  gain).  For  example,  amplitude  control  has  been 
demonstrated  through  carrier  injection*®’*''  or  depletion*®  in  a 
semiconductor  substrate,  both  of  which  shunt  the  capacitive 


region  of  the  metamaterial  elements.  However,  this  simple 
approach,  where  the  entire  capacitive  region  is  continuously 
covered  with  semiconductor,  is  insufficient  to  tune  the 
resonance  frequency. 

To  achieve  frequency  tunability,  due  consideration  must  be 
given  to  the  incorporation  of  the  semiconductor  into  the 
metamaterial  elements.  We  describe  hybrid  metamaterials  where 
precise  patterning  of  semiconductors  permits  frequency  tuning  of 
the  metamaterial  resonance  on  photoexcitation.  Our  results 
represent  an  important  initial  step  towards  extending  the  spectral 
range  over  which  a  specific  metamaterial  device  can  operate,  in 
addition  to  enabling  the  implementation  of  new  device  concepts. 

There  have  been  some  efforts  to  realize  frequency-tunable 
metamaterials*''"^*,  most  of  which  have  been  demonstrated  in  the 
microwave  frequency  region.  However,  it  should  be  noted  that 
these  designs  only  control  the  resonant  properties  of  one  or  a 
handful  of  individual  metamaterial  elements.  Thus,  these 
techniques  are  not  suitable  at  terahertz  (THz)  and  higher 
frequencies,  where  the  entire  metamaterial  composite  may  have 
more  than  1  x  10^  unit  cells  and  monolithic  integration  of  the 
tuning  elements  is  important. 

In  Fig.  la  we  show  a  scanning  electron  microscopy  (SEM) 
image  of  the  unit  cell  of  a  frequency- tunable  THz  metamaterial. 
It  is  a  variant  of  previously  demonstrated  designs®®"®'*  that  exhibit 
a  Lorentz-like  resonant  response  described  by  an  effective 
complex  permittivity  s{(o)  =  Sj  -|-  !E2,  where  Ej  and  E2  are  the 
real  and  imaginary  parts,  respectively.  The  split  gap  at  the  centre 
of  the  unit  cell  can  be  thought  of  as  a  capacitor,  where  charge 
accumulates  on  resonance.  By  constructing  the  capacitor  plates 
from  a  semiconductor  such  as  silicon,  we  can  control  the 
conductivity  in  response  to  the  photoexcitation  of  free  charge 
carriers,  thereby  altering  the  effective  size  of  the  capacitor  and 
tuning  the  capacitance,  C.  As  the  resonance  frequency  is  strongly 
dependent  on  the  capacitance,  that  is,  Wq  ~  (LC)"*/®,  where  L  is 
the  effective  inductance  of  the  split-ring  resonators  (SRRs),  then 
(Og  should  shift  monotonically  to  lower  frequencies  as  the 
photoexcitation  fluence  increases. 

The  metamaterials  were  fabricated  on  a  silicon-on-sapphire 
(SOS)  wafer.  The  (100)  silicon  layer  was  600  nm  thick,  with  an 
intrinsic  resistivity  >100  Cl  cm.  The  R-plane  sapphire  substrate  was 
530  p,m  thick.  The  planar  metamaterial  array  was  fabricated  using 
standard  photolithographic  methods,  consisting  of  electron-beam 
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Figure  1  Scanning  electron  microscopy  images  of  the  trequency-tunable 
planar  metamaterial,  a,  An  individual  unit  cell,  b,  A  periodically  patterned 
square  array.  All  dimensions  are  shown  in  micrometres  and  materials  are 
indicated  in  the  images.  The  polarization  of  the  normally  incident,  linearly 
polarized  THz  radiation  is  also  indicated  in  b.  H:  magnetic  field;  E:  electric  field. 


deposition  of  10  nm  of  titanium  and  200  nm  of  gold,  and  then  a  Uft- 
ofif  process.  The  metal  stripes  that  had  the  split  gaps  lay  along  the 
primary  flat  orientation  of  the  wafer,  that  is,  45°  counterclockwise 
from  the  projection  of  the  C-axis  on  the  R-plane.  Additional 
photoresist  was  patterned  by  photolithography  to  define  the 
silicon  regions  (capacitor  plates).  The  photoresist  patterns  and  the 
metallic  SRRs  together  served  as  a  mask  during  the  reactive  ion 
etching  (RIE)  of  the  unwanted  silicon  regions.  The  final  result 
(Fig.  la,b)  was  the  complete  removal  of  the  silicon  layer  except  for 
regions  under  the  metallic  SRRs  and  at  the  silicon  capacitor  plates. 
By  defining  the  SRRs  on  top  of  the  sihcon  layer,  we  avoided  any 
significant  steps  or  discontinuities  in  the  metallization. 

We  characterized  the  frequency-dependent  electric  response 
using  THz  time-domain  spectroscopy  (THz-TDS).  The  THz 
electric  field  was  coherently  measured  following  normal 
transmission  through  the  metamaterial  sample  or  a  bare  sapphire 
substrate  serving  as  the  reference.  The  THz  electric  field  was 
oriented  as  shown  in  Fig.  lb  to  drive  the  electric  resonance  of  the 
SRRs.  The  time-domain  data  were  temporally  windowed  to 
eliminate  the  effects  of  multiple  reflections  within  the  sapphire 
substrate.  Fourier  transformation  of  the  time-domain  waveforms 
then  provided  the  frequency-dependent  THz  electric-field 
amplitude  and  phase.  Dividing  the  sample  complex  spectrum  by 
the  reference  we  obtained  the  transmission  t{cj)  and  phase 
change  A(^(w)  through  the  metamaterial.  Photoexcitation  of 


Figure  2  Terahertz  electric-field  transmission  amplitude  of  tbe 
metamaterial  shown  in  Fig.  1.  a,  Experimental  measurements  as  a  function 
of  photoexcitation  power,  b,  Numerical  simulations  as  a  function  of 
silicon  conductivity. 


metamaterial  silicon  regions  was  achieved  using  optical-pump 
THz-probe  (OPTP)  spectroscopy^®.  The  OPTP  system  provided 
amplified,  ~30-fs,  near-infrared  laser  pulses  with  a  centre 
wavelength  of  800  nm  and  a  repetition  rate  of  1  kHz  to  optically 
excite  charge  carriers  across  the  1.12-eV  bandgap  of  silicon^®. 
The  photocarrier  lifetimes  were  orders  of  magnitude  longer  than 
the  picosecond  duration  of  the  THz  pulses.  This  enabled  a 
quasi-steady-state  characterization  of  the  metamaterials  as  a 
function  of  the  carrier  density,  which  could  be  adjusted  by 
changing  the  excitation  fluence.  In  our  measurements, 
photoexcitation  occurred  10  ps  before  the  arrival  of  the  THz 
pulses.  To  ensure  photoexcitation  uniformity,  the  pump  laser 
beam  was  expanded  to  a  diameter  of  ~1  cm  and  two  identical 
metal  apertures  of  3  mm  diameter  were  placed  in  front  of  the 
sample  and  reference.  The  focused  THz  beam  was  approximately 
3  mm  in  diameter  at  the  sample. 

The  experimental  results  for  the  frequency-tunable 
metamaterial  are  shown  in  Fig.  2a.  The  transmitted  THz  electric- 
field  amplitude  with  no  photoexcitation  is  shown  as  the  black 
curve.  A  transmission  minimum  t{a>)  =  19%  occurs  near 
resonance  at  1.06  THz,  whereas  off- resonance  transmission  is 
~90%.  As  the  photoexcitation  power  increases  (10-20  mW),  the 
resonance  initially  weakens  and  broadens,  showing  little 
perceptible  shift  in  frequency.  On  further  increase  in  the 
excitation  power  (50-100  mW),  the  resonance  continues  to 
weaken  and  saturates  at  t{(o)  =  49%,  but  now  shifts  to 
significantly  lower  frequencies.  Increasing  the  excitation  power 
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beyond  100  mW  causes  continued  shifting  to  lower  frequencies  and 
also  re-establishes  the  resonance  strength  and  narrows  the 
linewidth.  At  a  photoexcitation  power  of  500  mW  (a  fluence  of 
~0.5  mj  cm^^  or  carrier  density  of  ~2.5  x  10'®  cm^^)  we 
achieved  t(aj)  =  40%  at  850  GHz  (red  curve).  This  represents  a 
tuning  range  of  20%  in  the  resonance  frequency. 

To  explain  the  nature  of  the  frequency  tuning,  we  performed 
finite-element  simulations  using  commercial  software  Microwave 
Studios  by  CST  (ref  27).  We  simulated  a  single  cubic  unit  cell  as 
shown  in  Fig.  la,  with  the  SRR  located  at  the  centre.  The 
boundary  conditions  on  faces  parallel  to  the  propagation  vector 
k  were  set  as  perfect  electric  conductor  or  perfect  magnetic 
conductor  to  approximate  a  transverse  electromagnetic  plane 
wave  propagating  through  the  periodic  SRR  array  of  the  hybrid 
metamaterial  sample  or  the  reference  substrate,  as  shown  in 
Fig.  lb.  We  set  the  input  and  output  ports  on  the  top  and 
bottom  faces  of  the  cubic  unit  cell  in  the  air  and  sapphire, 
respectively,  which  minimizes  Fresnel  reflections  at  the  back 
interface  of  the  530-p,m-thick  sapphire  substrate.  In  simulations 
we  used  adaptive  mesh  refinement  to  ensure  an  accurate 
numerical  solution  with  a  short  simulation  time.  The  final  mesh 
typically  contained  about  50,000  tetrahedrons,  and  the  minimum 
edge  length  was  about  0.013  p,m.  The  simulation  produced  the 
frequency-dependent  complex  S-parameters  and  we  focus  here 
on  the  transmission  t{co)  =  |£,(w)/£;(w)|  =  |S2i(w)|  normalized 
by  a  blank  sapphire  substrate,  where  £,(w)  and  £;(w)  are  the 
transmitted  and  incident  fields,  respectively.  The  gold 
regions  were  simulated  as  lossy  metal  with  a  conductivity  of 
CTgoij  =  7  X  10®  Sm^',  and  the  silicon  underlying  the 

metamaterials  was  modelled  as  a  dielectric  with  85;=  11.7  and 
(Ts;  =  1  S  m^'.  The  sapphire  substrate  was  modelled  as  a  lossless 
dielectric  with  =  10.5,  the  average  value  of  e  along  the 

ordinary  and  extraordinary  axes^®.  The  photoexcited  silicon 
capacitor  plates  were  simulated  with  £5;=  11.7  and  a  pump- 
power-dependent  conductivity.  We  included  1-p.m  undercuts  in 
the  silicon  located  beneath  the  gold  regions,  an  attribute  of  our 
sample  caused  by  the  RIE  process.  Finally,  we  accounted  for  the 
irregular  geometry  (rounded  corners  and  chamfers  of  the 
fabricated  structure)  apparent  in  the  SEM  image  of  the  metamaterial. 

Simulation  results  are  shown  in  Fig.  2b,  and  excellent 
agreement  is  observed  both  quantitatively  and  qualitatively  when 
compared  to  the  experimental  measurements.  We  used  a  silicon 
conductivity  value  of  cr  =  50,000  S  m^'  to  simulate  a 
photoexcitation  power  of  500  mW  and  found  good  agreement 
with  experimental  data.  The  excited  charge-carrier  density  is 
proportional  to  the  excitation  power.  However,  the  conductivity 
increases  more  slowly  as  there  is  a  decrease  in  mobility  with 
photoexcited  carrier  density  due  to  increased  electron -electron 
scattering.  The  conductivity  values,  specified  in  Fig.  2b,  were 
found  to  be  in  reasonable  agreement  with  Drude  calculations  for 
photoexcited  silicon  with  the  fluences  specified  in  Fig.  2a. 

The  simulations  support  the  observed  experimental  behaviour. 
For  small  values  of  silicon  conductivity  there  is  an  initial  damping 
of  the  resonance  strength  due  to  ohmic  losses  in  the  silicon  plates, 
which  weakens  the  resonance  strength  as  shown  in  Fig.  2.  Further, 
the  small  silicon  conductivities  are  insufficient  to  enable  charge 
accumulation,  and  the  slight  change  in  dielectric  function  on  low 
fluence  photoexcitation  is  not  enough  to  alter  the  capacitance 
and  shift  the  resonance  frequency.  This  is  similar  to  the  results  of 
previous  studies*^’*''  in  which  switchable  THz  metamaterials  were 
created  by  incorporating  semiconductor  substrates.  Here  an  even 
lower  conductivity  was  required  to  completely  switch  off  the 
resonance  as  a  result  of  the  substrate  allowing  current  flow  across 
the  gaps  of  the  SRRs.  We  emphasize,  however,  that  in  the 
previous  studies'^"'®  the  injected  carriers  short-circuit  the  SRR 


□  ,  37  ^  □ 


Figure  3  Schematic  of  inductively  tuned  metamaterial  designs.  a,b,  The 
photoexcited  silicon  regions  form  parallel  current  paths  through  the  meandering 
loop  sections  (a)  and  across  the  regular  loop  sections  (b)  to  effectively  modify 
the  inductance  of  the  SRRs.  The  metal  and  silicon  regions  are  displayed  in 
orange  and  grey,  respectively.  The  dimensions  are  shown  in  micrometres. 


capacitance.  In  the  present  work,  the  patterned  semiconductor  is 
an  integral  part  of  the  metamaterial,  that  is,  the  SRR  capacitor 
plates.  As  the  conductivity  increases,  the  resonance  shifts  to  lower 
frequencies.  At  intermediate  conductivity  values,  more  charge 
accumulates  on  the  silicon  plates,  yielding  a  redshift  of  the 
resonance.  The  damping  remains  strong,  because  current  flow  is 
still  impeded  by  ohmic  losses  in  the  silicon.  Finally,  at  high 
conductivity  values,  the  silicon  begins  to  behave  more  like  a 
metal.  The  high  charge  carrier  density  imposes  little  resistance  to 
current  flow,  allowing  a  smooth  delivery  of  charge  between 
the  silicon  and  the  metallic  SRRs,  thereby  reducing  damping. 
The  charge  readily  accumulates  in  the  silicon,  contributing  to  the 
capacitance.  Effectively,  photoexcitation  changes  the  size  of  the 
capacitor  plates.  We  note  that  despite  the  reduction  in  oscillator 
strength  for  photoexcitation  powers  ranging  between  20  mW  and 
100  mW,  the  resonance  of  the  metamaterial  indeed  shifts  to  lower 
frequencies  as  expected,  from  1.06  THz  to  850  GHz,  a  tuning 
range  of  20%.  The  simulated  results  show  similar  behaviour,  with 
a  tuning  range  of  16%.  The  slight  discrepancies  in  frequency 
tunability  and  resonance  linewidth  between  the  experiment  and 
simulation  results  are  mainly  due  to  imperfections  in  sample 
fabrication,  which  were  not  fully  taken  into  account  in 
the  simulations. 

Eor  greater  flexibility  and  with  a  view  towards  applications  it 
would  be  ideal  to  have  additional  metamaterials  that  are  tunable 
to  higher  frequencies  on  photoexcitation.  In  Fig.  3  we  present 
simulated  studies  of  two  designs  that  accomplish  this  task  and 
yield  resonant  responses  in  the  THz  frequency  range.  Unlike  our 
previous  sample,  these  designs  use  frequency  tuning  by  altering 
the  inductance  L  of  the  metamaterial  element.  The  incorporation 
of  the  semiconductor  regions  shown  in  Fig.  3a  and  b  will,  under 
photoexcitation,  reduce  the  loop  length,  thereby  shifting  the 
resonance  to  higher  frequencies.  In  Fig.  4  we  show  the  simulated 
transmission  as  a  function  of  silicon  conductivity  for  both 
metamaterial  structures.  The  simulation  results  follow 
expectations:  the  resonances  shift  to  higher  frequencies  with 
increasing  semiconductor  conductivity.  As  with  our  previous 
design  in  Fig.  la,  the  resonance  initially  weakens  as  it  shifts,  but 
eventually  recovers.  Similar  tunabilities  are  achievable; 
simulations  show  (Oq  increasing  by  15%  in  Fig.  4a  and  40%  in 
Fig.  4b  for  the  structures  in  Fig.  3a  and  b,  respectively.  However, 
it  should  be  noted  that  the  conductivity  values  needed  to  realize 
such  tunabilities  are  one  order  of  magnitude  higher  than  those 
used  in  Fig.  2b.  This  results  from  the  significant  amount  of 
current  required  to  traverse  the  silicon  regions  to  reduce  the 
effective  loop  size.  Thus  the  silicon  conductivity  has  to  be 
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Frequency  (THz) 


Figure  4  Simulated  THz  electric-field  transmission  amplitude  of 
metamaterial  designs  presented  in  Fig.  3.  a,b,  Simulations  of  the  structure  of 
Fig.  3a  (a)  and  the  structure  of  Fig.  3b  (b),  both  using  the  silicon  conductivity 
values  shown  in  the  key. 


comparable  to  that  of  gold,  which  is  difficult  to  achieve  by 
photoexcitation  of  silicon.  In  contrast,  there  is  much  less  current 
travelling  through  the  silicon  regions  used  to  construct  the 
capacitor  plates  in  Fig.  1.  In  any  case,  the  simulated  structures 
along  with  the  experimentally  demonstrated  frequency-agile 
metamaterials  (Figs  1  and  2)  will  be  useful  for  the  construction 
of  future  novel  devices.  The  extended  functionality  offered  by 
the  hybrid  metamaterials  presented  here  will  undoubtedly  play 
a  key  role. 
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